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ABSTRACT 

Recent discoveries of compact (sizes < R©) debris disks around more than a dozen of metal-rich white 
dwarfs (WDs) suggest that pollution of these stars with metals may be caused by accretion of high-Z 
material from the disk. But the mechanism responsible for efficient transfer of mass from a particulate 
disk to the WD atmosphere has not yet been identified. Here we demonstrate that radiation of the 
WD can effectively drive accretion of matter through the disk towards the sublimation radius (located 
at several tens of WD radii), where particles evaporate, feeding a disk of metal gas accreting onto the 
WD. We show that, contrary to some previous claims, Poynting-Robertson (PR) drag on the debris 
disk is effective at providing metal accretion rate Mpji ~ 10^ g s~^ and higher, scaling quadratically 
with WD effective temperature. We compare our results with observations and show that, as expected, 
no WD hosting a particulate debris disk shows evidence of metal accretion rate below that produced 
by the PR drag. Existence of WDs accreting metals at rates significantly higher than Mpn suggests 
that another mechanism in addition to PR drag drives accretion of high-Z elements in these systems. 

Subject headings: White dwarfs — Accretion, accretion disks — Protoplanetary disks 



1. INTRODUCTION. 

Recent infrared observations with Spitzer and ground- 
based facilities (Zuckerman & Becklin 1987; Graham et 
al. 1990; Farihi et al. 2010) revealed near-infrared ex- 
cesses around more than a dozen of metal-rich white 
dwarfs (WDs). This emission was interpreted as repro- 
cessing of the WD radiation by refractory material re- 
siding in an extended, optically thick and geometrically 
thin disk (Jura 2003; Jura et al. 2007), similar to the 
rings of Saturn (Cuzzi et al. 2010). 

It was hypothesized by Jura (2003) that such compact 
disks of high-Z material may be naturally produced by 
tidal disruption of asteroid-like bodies entering the Roche 
radius of the WD. This idea naturally explains the well 
defined outer radii Rout ~ 1 R© of disks since the Roche 
radius Rp for tidal disruption of a self-gravitating object 
of normal density (p 1 g sm~^) by Mi, ^ central 
mass is Rr - {M^^lpf/^ w 1 R© (Jura 2003). 

Moreover, Jura (2003) went on to suggest that the 
high-Z material contained in compact disks is respon- 
sible for the observed metal enrichment of a significant 
fraction of WDs. Accretion of this circumstellar mate- 
rial at high enough rate Mz onto WD could maintain 
a non-zero abundance of metals in the WD atmosphere 
against rapid gravitational settling. This scenario thus 
provides a promising alternative to the previously widely 
discussed interstellar a,ccvetion model of WD metal pollu- 
tion (Dupuis et al. 1993), which is known to have serious 
problems. 

Observed abundances of heavy elements in WD at- 
mospheres and theoretical calculations of their gravita- 
tional settling imply typical metal accretion rates Mz ~ 
_ iQio g (Farihi et al. 2009, 2010). An evolving 
disk of debris must be able to supply such high Mz to 
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the WD. However, the question of how high-Z elements 
get transported to the WD atmosphere from the ring of 
solid particles, which does not extend all the way to the 
WD surface, has not yet been answered. 

In principle, a dense ring of particles should evolve 
simply because of the angular momentum transfer due 
to inter-particle collisions, in full analogy with the rings 
of Saturn. However, the evolution time scale of Saturn's 
rings due to this process is too long, ~ 10^ yr (Salmon 
et al. 2010), and resultant values of Mz are negligible. 

Another natural mechanism driving debris towards the 
WD is due to stellar radiation interacting with the disk 
and giving rise to the Pointing-Robertson (PR) drag 
(Burns et al. 1979). Previously, Farihi et al. (2010) 
claimed that PR drag cannot provide Mz higher than 
10^ — lO'* g s^^ and dismissed this process as irrelevant. 
The goal of this paper is to critically re-examine the ef- 
fect of radiative forces on the debris disk evolution, and 
to show in particular that PR drag can give rise to Mz 
inferred from observations. 

2. MASS ACCRETION DUE TO RADIATIVE FORCES. 

We envisage the following conceptual picture of the 
circum-WD environment. A dense disk (or ring) of par- 
ticles lies inside the Roche radius Rr of the WD and 
evolves under the action of external agents, e.g. radia- 
tion forces. Particles migrate through the disk towards 
the sublimation radius i?^, where their equilibrium tem- 
perature equals the sublimation temperature T^: 
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where i?^ is the WD radius, Tg « 1500 K for silicate 
grains, and r^,4 = T^/IO^ K is the normalized stellar 
temperature T^. Taking R^ « O.Oli?© typical for mas- 
sive (Af^ > 0.6 Mq) WDs (Ehrenreich et al. 2011), one 
finds Rs ~ 0.2 Kq, in agreement with observationally 
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inferred inner radii of compact debris disks (Jura et al. 
2007, 2009a). Thus, for cool WDs Rout/Rin is several. 

Particles sublimate at Rs feeding a disk of metallic gas, 
which we assume to be transparent to stellar radiation 
(gaseous component has been detected in several WDs 
with debris disks, see Gansicke etal 2006). Gaseous disk 
viscously evolves, extending all the way to the WD sur- 
face and providing means of metal delivery from Rs to 
the star. Assuming conventional a-parametrization of 
viscosity (Shakura & Sunyaev 1973) v = ac^/flK, the 
characteristic viscous time in the disk of metallic gas is 
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where M^, l = M^^/Mq, and /i28 is the mean molecular 
weight of the metallic gas normalized by 28mp (value of 
/i for pure Si). This timescale is short enough for Mz to 
be determined predominantly by evolution of the disk of 
solids. 

We take the disk of particles to be geometrically thin 
and characterize it at each point by optical depth r de- 
fined as 

_ 3 Ed 
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where is the surface mass density of debris disk, and p 
and a are the bulk density and characteristic size of the 
constituent particles (we leave a unspecified for now). 
Observations suggest (Jura 2003; Jura et al. 2007) that 
T > 1 but we leave t arbitrary. 

Particles in a narrow annulus just outside the sublima- 
tion radius are exposed to direct starlight at normal inci- 
dence and are heated to sublimation temperature. Farihi 
et al. (2010) studied the effect of PR drag on particles in 
this region. Here we first explore the effect of radiative 
forces on the rest of the disk, which has only its surface 
illuminated by the central star at grazing incidence. To 
account for the finite size of the WD we adopt a "lamp- 
post" model of stellar illumination, in which radiation is 
emitted by two point sources of luminosity L*/2 each lo- 
cated at height Z = {4:/3^T)R^, above and below the disk 
plane, see Figure [T] Then 
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is the incidence angle of radiation on the disk surface at 
distance r from the star, and a 1 since r > i?^ 3> i?^. 
The one-sided incident energy flux per unit surface area 
of the disk is Fr = (L^,/67r^)(i?^,/r^), in agreement with 
Friedjung (1985). 

Radiation falling on the disk gives rise to two forces 
that can change angular momentum of the disk: PR 
drag and the Yarkovsky force. The latter owes its ex- 
istence to thermal inertia of spinning objects that causes 
re-emission of absorbed energy in the direction different 
from that of incoming radiation (Bottke et al. 2006). Az- 
imuthal components of these forces per unit area of the 
disk can be represented as 
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where c is the speed of light, and V' is the factor that 
characterizes the phase lag of radiative momentum de- 
position in disk particles and is specified later in fj223] 




Disk of solids 



Fig. 1. — "Lamp-post" model of debris disk illumination by the 
WD at small incidence angle a. WD emission is represented by 
two point sources at height ±Z to account for the finite size of the 
WD. An inner edge of the disk of solids (black) near the sublimation 
point at Rs is shown, including an annulus of directly exposed solid 
material (gray). White arrows indicate accretion of solids, black 
arrows indicate accretion of metal gas inside of . Disk thickness 
h <C R* is not shown to scale. Inset on top schematically shows 
the run of the disk optical depth with radius. See text for details. 



and §2 2.2l for PR drag and Yarkovsky force respectively. 
Factor (/>r characterizes the efficiency of radiative momen- 
tum absorption by the disk surface. In the geometrical 
optics limits (particle size a much larger than character- 
istic wave length of stellar radiation) 

(/)^ « 1 - e^^ii , T||=a"V, (6) 

where T|| is the optical depth encountered by incident 
photons as they traverse the full disk thickness. Optically 
thick debris disks have (pr — ^■ 

Azimuthal force causes radial drift of disk material at 
speed Vr = 2/^/(r2i<-I]c;), where ^Ik is the Keplerian an- 
gular frequency. This gives rise to mass transport at the 
rate 



(7) 



Note that for radiative forces both f^p and Mz are inde- 
pendent of Erf. 



2.1. Poynting- Robertson drag. 
In the case of PR drag factor ip is (Burns et al. 1979) 
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where we used equations ^ and ([7]). 

Numerical estimate in equation ^ does not agree with 
the calculation of Farihi et al. (2010) and we demonstrate 
the origin of this discrepancy by following these authors 
and considering the inner edge of the disk where particles 
sublimate (at r = i?^). At the edge there is an annulus 
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of the disk directly exposed to starlight (i.e. T|| _(j < 1, see 
Figure [l]) with radial width 
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which follows from the relationship Ta/h ^ ''"||,a/^ (here 
h <^ Ri, is the disk thickness, and Ty < 1 are the 
vertical and horizontal optical depths of the annulus). 
Particles inside this annulus experience PR drag as in- 
dividual objects, with little shadowing from nearby par- 
ticles. The characteristic decay time of their semimajor 
axes is (Burns et al. 1979) 
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However, it only takes the fully exposed particles time 
tcross ^ tpR{l/ Rs) to cross the annulus, reach the sub- 
limation point, and evaporate, at which point they stop 
shadowing objects behind them. Dividing the mass 



2nRslT,d contained within the annulus by tc 
the mass accretion rate inside the annulus: 
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where we used definition ([3]). Mass conservation requires 
that Ma = MpR,s = Mpb{Rs) so that the optical depth 
inside the annulus is Ta ^ as4>r ^ ^ {as = oi{Rs)) and 
from equation (|10|) the width of the annulus is h/ {as<j)r)- 
Thus, as long as the disk outside the annulus is not trans- 
parent to incident stellar radiation (i.e. T|| > 1 or r > a) 
T drops from its value Tg outside of i?s -I- / to inside. 
This is because fully exposed particles in the annulus ex- 
perience stronger PR drag and migrate towards the WD 
faster than particles outside the annulus. The transi- 
tion occurs at the distance / outside the disk edge, where 
T\ia ~ 1, see Figure m 

From observational point of view, the most interesting 
characteristic of the disk is the value of MpR at subli- 
mation point i?s, since this is the rate at which disk of 
solids feeds gaseous disk, from which mass is transported 
onto the WD atmosphere by viscous torques. Evaluating 
expression © at r = i?s using equation ([1} we find 
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where i?^_2 = R*/10^^ R©; then Mpp{r) = 
Mpp^siRs/r). Note that Mpp^g depends on only 
quadratically. 

Previously, Graham et al. (1990) obtained an empirical 
estimate of Mpp by assuming accretion to be driven by 
PR drag and dividing the observationally inferred disk 
mass (surface area multiplied by pa, taking disk to be a 
monolayer of particles) by tpp, see equation ()11|) . Even 
though their numerical estimate of Mpp « 1.8 x 10* 
g is not very different from reality this method of 
deriving MpR is not well motivated. Later, Farihi et 
al. (2010) found a value of Mpp < 10^ g s""'^ much lower 
than MpR^s predicted by equation because they used 
tpp instead of tcross ^ tpp in computing Ma- 



2.2. Yarkovsky force. 

If circum-WD debris disks behave similarly to dense 
planetary rings one expects collisions in the disk to align 
particle spins predominantly parallel or anti-parallel to 
the disk normal, with non-zero and positive average spin 
(prograde mean rotation) and significant random spin 
component (Salo 1987; Ohtsuki & Toyama 2005). The 
characteristic spin frequency is cj ^^k- Thermal inertia 
of spinning particles causes re-emission of absorbed stel- 
lar energy in the direction different from the radial. This 
gives rise to azimuthal Yarkovsky force, which causes out- 
ward orbital migration of disk particles spinning in pro- 
grade sense (Bottke et al. 2006). We now check whether 
this force can affect the PR-driven inward accretion of 
solids by estimating its typical magnitude. 

Yarkovsky phase lag factor ipy depends on two dimen- 
sionless parameters. One of them R' = a /In is the ratio 
of particle size a to the penetration depth of the ther- 
mal wave li, — K/ pCpijj, where K and Cp are the 
thermal conductivity and specific heat of particle ma- 
terial. Another is the thermal parameter (roughly the 
ratio of thermal time at depth ly to object's spin period) 
6 = ^ KpCpLo I gT'^ , where T is the surface temperature 
of the body. The actual dependence of V'y on R! and 
8 has been calculated in Vokrouhlicky (1998, 1999) and 
these references should be consulted for details. For our 
current purposes it suffices to note that phase lag factor 
takes on a maximum value i/'r, max 

0.1 when R' > 1 

and ~ 1 , and is smaller for other values of R' and 8 . 

To evaluate l^, and 8 we assume that disk particles 
have composition typical for terrestrial silicate minerals 
such as olivine, in agreement with elemental abundances 
of metal-rich WD atmospheres (Zuckerman et al. 2007; 
Klein et al. 2010). Then if « 1.5 x 10^ erg s^^ cir-^ 
K-i and Cp w 10^ erg g^^ K'^ at T > 10^ K (Roy et 
al. 1993). Taking uj = fix we can write 

where l^^g and 8^ are the values of li, and 8 at sublima- 
tion radius r — Rs'. 
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Particle size a needed for calculation of R' is not con- 
strained by observations, although Graham et al. (1990) 
suggested that debris disk may consist of 10 — 100 cm 
particles. If a is several cm or larger (so that R' > 1) 
then ipY,max ^0.1 for fully illuminated particles at Rs, 
see Fig. 3 of Vokrouhhcky (1998). 

However, in the dense debris disk only the upper (or 
lower) parts of particles, close to their spin axes, are di- 
rectly illuminated by anisotropic starlight. This is be- 
cause the latter illuminates the disk at small angle a and 
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most of particle surface is shadowed by other particles 
nearby (thermal emission of these particles is on average 
isotropic in horizontal direction). From heuristic argu- 
ments one expects stellar heating to be mainly deposited 
at small colatitude ^ a^/^ ^ 1 from the particle spin 
axis. It is easier for thermal conduction to isotropize sur- 
face temperature distribution around the spin axis over 
this small polar cap region rather than over the whole 
particle surface (as is assumed in standard calculation of 
4'y), and this additionally lowers ipy- 

Thus, even if material properties of disk particles are 
favorable for maximizing ipY (as is the case for our es- 
timates here) one still should expect ipy ^ 0.01. Equa- 
tions (U) and (O result in the following expression for 
Yarkovsky-induced outward mass accretion rate at i?^: 
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This is somewhat higher than Mpn^s, thanks to the high 
adopted ^py = 10^^ exceeding ippR, see ([5]). However, 
this estimate may be very optimistic because of our poor 
knowledge of material properties and sizes of disk parti- 
cles. 

Finally, we note that observations do not reveal the 
existence of metal-poor WDs with debris disks around 
them, which could plausibly exist if outward migration 
of particles due to the Yarkovsky force were effective at 
stopping the PR-driven inward mass accretion. This pro- 
vides additional evidence that Yarkovsky effect does not 
play significant role in the circum-WD debris disk evolu- 
tion. 

3. COMPARISON WITH OBSERVATIONS. 

We now compare our theoretical predictions with data 
on Mz inferred from observations'^_ of metal-rich WDs. 
In Figure [2] we plot the values of Mz from Farihi et al. 
(2009, 2010) vs. stellar effective temperature T,, sepa- 
rately for WDs with and without debris disks detected 
via IR excesses. We also plot our analytical prediction 
for Mpfi^s{Ti,) for different values of i?* and Tg. 

It is easy to see from Figure [2] that all systems with de- 
tected debris disks exhibit Mz higher than Mpus- This 
is very encouraging since whenever WD is orbited by par- 
ticulate disk one expects PR drag to set a lower limit of 
the mass flux at the level of Mpp s- This lower bound on 
Mz is hard (if at all possible) to avoid. This constraint 
naturally explains a strong positive correlation between 
Mz and the rate of disk occurrence found by Farihi et al. 
(2010): the fraction of stars with compact debris disks 
is significantly higher for WDs with high Mz simply be- 
cause systems with disks cannot have low Mz < Mpp^s- 
All this strongly suggests that the PR drag indeed plays 
important and visible role in circum-WD debris disk evo- 
lution. 

It is worth emphasizing that the lack of disk-hosting 
systems at low Mz is not caused by some observational 

^ These estimates are somewhat model-dependent, see Jura et 
al. (2009b). 
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Fig. 2. — Mass accretion rates of high-Z elements Mz inferred 
from elemental abundances measured in WD atmospheres. Metal- 
rich WDs without debris disks around them (with no detected IR 
excess) are depicted as crosses, WDs with disks are displayed by 

open points. Curves show analytical prediction for MpE.,s (for 
different WD radii Ri, and sublimation temperatures Ts ) given by 
equation I I13I I. Note that all WDs with disks tend to lie above 
^'IpR,s curve, implying that whenever a debris disk is present the 
Pointing-Robertson drag sets a lower limit on the metal accretion 
rate. 



bias: disk-bearing WDs with Mz ^ Mpp,s such as GD 
16 (Farihi et al. 2009), GD 56 (Jura et al. 2009a), G166- 
58 (Farihi et al. 2008) have (sometimes highly) significant 
detections of IR excesses by Spitzer. 

It is also clear from Figure [2] that there are metal- 
rich systems both with and without disks, in which Mz 
significantly exceeds Mpp,s for a given T*, sometimes by 
several orders of magnitude. The most dramatic case is 
that of DAZB WD GD 362 which has Mz « 2.5 x lO^o g 
s^^ (Farihi et al. 2009), exceeding Mpp^s corresponding 
to its = 10, 500 K by a factor of > 300. Existence of 
such objects clearly implies that at least in some cases 
an additional accretion mechanism must operate on top 
of the PR drag, giving rise to very high Mz- We suggest 
such a mechanism that naturally operates in presence of 
a debris disk in Rafikov (2011, in preparation). 

Systems without detected debris disks tend to occupy a 
broad range of Mz, both above and below Mpp^s- They 
may be interpreted as WDs that were orbited by com- 
pact debris disks in recent past but have now completely 
(or largely) lost their disks to accretion. Without active 
external source metals sediment out from atmospheres of 
these WD resulting in lower inferred Mz in systems that 
had longer time since their debris disk disappearance. 

Lifetime of a debris disk of mass Md accreting via PR 
drag alone can be estimated as M^/Mpp « 2 — 4 Myr 
for Md = 10^2 g (roughly corresponding to the mass 
of 200 km asteroid) and Mpp typical for a WD with 
Ti, = 10"* K. This estimate only marginally agrees with 
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observational constraints suggesting several 10^ yr disk 
lifetime (Farihi et al. 2009). However, one should keep 
in mind that this is an upper limit on the lifetime that 
is realized if no other processes apart from the PR drag 
drive accretion of solids through the disk. 

4. DISCUSSION. 

In conclusion we would like to emphasize the robust na- 
ture of the mass accretion due to PR drag. Our estimate 
(fT3|l of Mpji^s depends neither on the disk properties 
such as Srf (as long as r|| > 1) or nor on the mate- 
rial properties of its constituent particles. Also, MpB.,s 
does not vary much as or R^, change within reasonable 
limits. 

One may worry that our calculation essentially dis- 



regards the details of the particle sublimation process. 
However, Mpp^s is set by stellar illumination outside the 
innermost annulus of directly exposed particles at the 
inner edge of the disk, see H2 2.11 There sublimation is 
irrelevant and our calculations are robust. Interior of 
that point Mz does not change by continuity and, as a 
result, mass accretion rate ends up being insensitive to 
how exactly directly exposed particles sublimate. 
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